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Abstract: Forthcoming projects such as the Dark Energy Survey, Joint Dark En- 
ergy Mission, and the Large Synoptic Survey Telescope, aim to measure weak lensing 
shear correlations with unprecedented accuracy. Weak lensing observables are sensi- 
tive to both the distance-redshift relation and the growth of structure in the Universe. 
If the cause of accelerated cosmic expansion is dark energy within general relativity, 
both cosmic distances and structure growth are governed by the properties of dark 
energy. Consequently, one may use lensing to check for this consistency and test 
general relativity. After reviewing the phenomenology of such tests, we address a 
major challenge to such a program. The evolution of the baryonic component of 
the Universe is highly uncertain and can influence lensing observables, manifesting 
as modified structure growth for a fixed cosmic distance scale. Using two proposed 
methods, we show that one could be led to reject the null hypothesis of general 
relativity when it is the true theory if this uncertainty in baryonic processes is ne- 
glected. Recent simulations suggest that we can correct for baryonic effects using a 
parameterized model in which the halo mass-concentration relation is modified. The 
correction suffices to render biases small compared to statistical uncertainties. We 
study the ability of future weak lensing surveys to constrain the internal structures of 
halos and test the null hypothesis of general relativity simultaneously. Compared to 
alternative methods which null information from small-scales to mitigate sensitivity 
to baryonic physics, this internal calibration program should provide limits on devi- 
ations from general relativity that are several times more constraining. Specifically, 
we find that limits on general relativity in the case of internal calibration are de- 
graded by only ~ 30% or less compared to the case of perfect knowledge of nonlinear 
structure. 
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1. Introduction 

Current data indicate that the cosmological expansion is accelerating (e.g., Refs. ^ 
1^, H, I, 1^, D, 1^, H, ^). This accelerating expansion is one of the most profound 
discoveries of the past decade and points to a fundamental gap in our understanding 
of the universe. The cause of the accelerating expansion has been dubbed dark 
energy. Indeed dark energy is simply a name for the fundamental puzzle that is the 
nature of cosmic acceleration and an enormous amount of effort is now being devoted 
to shedding light on the causative agent of cosmic acceleration. In broadest terms, 
two options have been put forth to describe accelerated cosmic expansion. The first 
is that cosmic acceleration is caused by some as yet unidentified contribution to the 
stress energy of the universe. This option includes vacuum energy (observationally 
indistinguishable from a cosmological constant) and dynamical models of dark energy. 
The second option is that gravity deviates from the general relativistic description 
on large scales, an option we refer to as modified gravity. 

Viable options to general relativity for which definite predictions have been made 
are few and far between, so many authors have suggested that a fruitful way to apply 
forthcoming data will be to check for the mutual consistency of different observable 
phenomena with the predictions of general relativity jT^, [TT|, |T^, |T^, |T^]. 
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The basic idea is that one may obtain observational handles on both the distance- 
redshift relation and the growth rate of cosmic expansion. Within general relativity, 
both of these can be predicted from the same dark energy parameters, so it is possible 
to measure, for example, distance and then check for consistency with the growth 
of cosmic structure. Forthcoming weak gravitational lensing surveys will provide 
the most powerful means to probe the matter distribution in the universe directly 
and are an indispensable piece of any such consistency check [TU, T^\. A significant 
literature already exists detailing the power of tomographic weak lensing to constrain 
dark energy parameters, including Refs. [jl8|, O, EO, |22|, I 



24, m 



27, 28, 29 



30| , pl[] . In fact, the effectiveness of any given matter fluctuation to serve as a lens 
is sensitive to the distance scale of the universe and so weak gravitational lensing 
alone will provide a very powerful check for the consistency of general relativity. 
In this manuscript, we study the effectiveness of such a consistency check in light 
of recent studies that indicate that predictions for weak lensing observables will 
contain significant inherent uncertainties due to the poorly-understood behavior of 
the baryonic component of the universe p2|, p3|, 531 EB, pll]. 



The strategy for checking the consistency of general relativity and searching for 
a sign of modified gravity outlined in the previous paragraph seems simple enough. 
However, making predictions for weak lensing power spectra is fraught with numer- 
ous practical difficulties. Weak lensing will use information on scales where density 
fluctuations are well beyond the linear regime (relative overdensities 5 > 1). To 
utilize this information at the level that forthcoming observational programs will 
permit, numerical simulations must be able to predict the nonlinear matter power 
spectra to better than a percent ^ . Sufficient precision should be achievable by 
brute force using dissipationless, A^-body numerical simulations that treat only the 
evolution of the cosmic density field under gravity only (e.g., Refs. 



38, 



The most challenging obstacle to such precise predictions is the significant infiuence 
that baryonic processes have on matter power spectra on scales of interest. This has 
been pointed out in both analytic 



J2|, ^ and numerical studies [0, p5|] . 
The problem posed by baryonic processes is severe because they cannot be mod- 
eled directly and all such calculations rely on relatively poorly-constrained, effec- 
tive models that approximate the net, large-scale, effects of processes that occur on 
scales far below the numerical resolution that may be achieved with any simula- 
tion. Ref. [BH] shows that the influences of baryons may be large compared with the 
statistical uncertainties expected of future surveys and that different treatments of 
baryonic physics lead to notably different matter power spectra on relevant scales. 
Ref. |31[| studied a parameterized model for baryonic influences that could treat all 
of the simulations of Ref. and showed that such a parameterized model could 
be calibrated self-consistently using weak lensing data alone, yielding meaningful 
constraints on both the dark energy and the effective models of baryonic physics. 
However, extending the results of Ref. PT| to consistency checks of general relativity 
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requires care. Baryonic processes manifest as an inability to predict the evolution 
of the density field even for a fixed cosmic expansion history. This is, in part, what 
Ref. 1^ relied upon in their self-calibration program and this is precisely what pro- 
posed tests of the consistency of general relativity rely upon. 

The strategy behind all efforts to study dark energy using gravitational lensing 
relies on the ability to produce reliable A^-body simulations. In fact, most stud- 
ies implicitly assume such a simulation campaign will be performed prior to any 
data analysis and proceed to estimate the power of forthcoming experiments using 
contemporary fitting formulae. In fact these fitting formulae are not sufficiently 
precise to treat forthcoming data, they may be subject to fundamental limitations 



1^, ^ and, aside from a handful of specific cases, they have not been 

generalized to treat dark energy. A numerical campaign will be necessary to address 
dark energy and addressing observables using dissipationless A^-body simulations is 
a challenging, but tractable problem [^, ^ . The results of A^-body simulation cam- 
paigns represent precise solutions to the idealized problem of computing the density 
or lensing fields in a variety of cosmologies absent baryonic physics. As with dark 
energy constraints, the strategy of self-calibrating the net influences of baryons is 
based on the assumption that such a set of simulations without baryons are available 
and that a simple set of prescribed corrections can be applied to the A^-body results 
to describe baryonic effects. The reason is that direct numerical calculations that 
treat the physics of baryons are not achievable in the foreseeable future, given current 
computational limitations. This is the context of the present paper. 

In this paper, we extend the results of Ref. |31[| to address consistency checks of 
general relativity. We present a brief pedagogical discussion of the manner in which 
independent modifications to the distance scale of the universe and the growth rate of 
structure produce non-degenerate changes in weak lensing statistics, enabling weak 
lensing alone to serve as a powerful consistency check for general relativity. We 
then move on to discuss the interplay between modifications to convergence spectra 
caused by baryonic processes and the dark energy. We show explicitly that our cur- 
rent ignorance regarding the evolution of the baryonic component of the universe, 
its influence on the evolution of inhomogeneities of dark matter, and the process of 
galaxy formation places a severe limitation on our ability to test the consistency of 
general relativity. To be specific, we consider two tests. In one test, the growth of 
structure and the cosmic distance scale are assumed to arise from two different effec- 
tive equations of state Wg and Wd, and we test our ability to rule out the hypothesis 
that the two are equal as they would be in general relativity. If unaccounted for, 
our limited ability to predict convergence spectra can lead to biases that drive these 
two parameters to disagree at levels as large as ~ 8a (depending upon details) when 
general relativity is the correct description of gravity. We perform a similar test on 
the gravitational growth index parameter 7, introduced in Ref. and likewise find 
significant biases if baryonic effects are not treated. We show that it is possible to 
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eliminate these biases by disregarding the small-scale shear information from forth- 
coming surveys, but the cost is a factor of ~ 2 — 4 degradation in the constraining 
power of such surveys. Lastly, we study the ability of forthcoming surveys to test 
general relativity and constrain baryonic processes simultaneously within a single 
data set. We show that this is a promising option as the biases can be eliminated at 
a cost of only ~ 20% in parameter constraints. This implies that significant effort 
should be devoted to developing robust parameterizations of the nonlinear evolution 
of structure as we begin to realize the next generation of imaging surveys. These 
methods should extend current techniques by including realistic descriptions of the 
effects of baryons that are not so general as to be completely arbitrary but do allow 
sufficient freedom so as to reflect our ignorance of the influence of the baryonic sector. 

In the following section, we describe our methods, including our parameteriza- 
tions of cosmological expansion and structure growth. In § |], we present the results 
of our study. We begin in § |3]l| by illustrating the power of these consistency checks. 
This section contains a compilation of facts and figures found dispersed throughout 
the existing literature. In § p.2| , we show how neglect of baryonic processes may lead 
to biases that would mimic an inconsistency in general relativity in simpler treat- 
ments of weak lensing observables. In § |3.3| , we present our results that show that 
additional parameter freedom may be added to predictions of the growth of structure 
to account for the influence of baryons and make such consistency checks robust to 
nonlinear processes. In § we discuss our results and we summarize our efforts in 
§1- 



2. Methods 



2.1 Weak Lensing Observables 



We treat weak lensing observations as a consistency check of general relativity. Aside 
from our parameterization of dark energy, we perform all calculations as described in 
detail in Ref . |^ . We summarize those calculations here and refer the reader to this 



reference for more detail. We consider number density-weighted convergence power 
spectra and cross-spectra from A^tom = 5 tomographic redshift bins as our primary 
observables: 

'H{z)Dl{z) 

The five photometric redshift bins are spaced equally in redshift between a minimum 
photometric redshift of 2;™'° = and a maximum z^^^ = 3. Previous work showed 
that parameter constraints are saturated with iVxoM = 5 and that further binning is 



^^(^) = / J'.'Ji.i Pik = i/Dj,, z). (2.1) 



unnecessary [^, We have verified that this remains so in the models we consider. 
In Eq. ( |2.1| ), is the angular diameter distance, H{z) is the Hubble expansion rate, 
P{k,z) is the matter power spectrum at wavenumber k and redshift z, and 
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are the lensing weight functions. Given the true redshift distribution of sources in 
bin i, drii/dz, the lensing weight is 

(2.2, 

We assume that the true redshift distribution of sources is 

dn{z)/dz = Aniz/zof exp[-iz/ zof]/ ^2^„ (2.3) 

where n represents the total density of source galaxies per unit solid angle, which 
varies from survey to survey, and we assume Zq ~ 0.92 so that each survey has 
a median redshift of z^ed = 1- We assume that the probability of a photometric 
redshift Zp given a spectroscopic redshift z is Gaussian with a mean value of Zp given 
by z (no bias) and a dispersion = 0.05(1 + z), and compute drii/dz for each 
photometric redshift bin as described in Ref. |45[] . 

We estimate the constraining power of forthcoming weak lensing surveys using 
the Fisher information matrix. Indexing the observables of Eq. (|2.1| ) by a single label, 
we write Oa = P^^, where each i,j map onto a unique A, and the Fisher matrix of 
weak lensing observables is 

= 'f,(2e + 1)/.,, J: ^|C-1.b^ + Fj,. (2.4) 

^min A,B 

The Pa represent the parameters of the model, Greek indices label model parameters, 
capital Latin indices label unique observables, and lower-case Latin indices label 
tomographic redshift bins, /sky is the fraction of sky covered by an experiment. Cab is 
the covariance matrix of observables, and the sum begins at imm = 2/ \/ fsky and runs 
to some ^max- We generally take £max = 3000 to ensure that we only use scales where 
assumptions of weak lensing and Gaussian statistics are valid ^ ^ pO| 



However, we present many of our primary results (Figures ^ |], 0) as a function 
of ^max so results for alternative choices are simple to extract. Observed spectra 
contain both signal and noise, = P^^ +ni(5ij(7^), where ni is the surface density of 
source galaxies in bin i and (7^) is the intrinsic source galaxy shape noise. We follow 
convention in setting (7^) = 0.2 and allowing differences in shape noise between 
different observations to be absorbed into Ui. The covariance between observables 
P'j and P^^ is Cab = P^P^ + P^P^k-: where i and j map to A, and k and 1 map to B. 

The inverse of the Fisher matrix is an estimate of the parameter covariance near 
the maximum of the likelihood. The measurement error on parameter a marginalized 
over all other parameters is 

a(p„) = [F-Xa. (2.5) 

The second term in Eq. ( |2.4| ) incorporates Gaussian priors on model parameters. 
We assume modest priors on several cosmological parameters individually, so that 
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F^p = Sais/cTa^ wheie aa is the assumed prior uncertainty on parameter a. Un- 
less otherwise stated, when we refer to the uncertainty in a parameter or subset of 
parameters, we are referring to the uncertainty in the parameters under discussion 
after marginahzing over the remaining parameters of the model. This formalism also 
provides an approximation for biases in cosmological parameter estimators due to 
unknown, untreated, systematic offsets in observables. Let AOa be the difference 
between the true observable and the prediction for that observable absent the sys- 
tematic. The induced bias in the estimator of parameter a due to the neglect of the 
systematic offset is 

do 

6p^ = 5^[F-i]„^5^(2£+ l)/3ky5^ AOA[C7-i]ABy^. (2.6) 

p e A,B 

We model the matter power spectrum using the phenomenological halo model 
0, |5^ (for a review see Ref. particular implementation of 



Ref. ^T|. This model utilizes standard fitting formulae for halo abundance and halo 



bias 1 5^. It is known that existing fitting formulae are not yet sufficiently precise 



to treat forthcoming data sets but it is likely that this can be overcome. 

Our approach is premised on the idea that dissipationless simulation programs will be 
carried out to calibrate these quantities to the necessary precision (or some equivalent 
strategy that utilized the simulation data directly), and that baryonic processes are 
the only processes that are not treated with sufficient precision. The non-standard 
modification to the halo model that we consider concerns the distribution of matter 
within dark matter halos. Typically, it is assumed that on average the mass density 
within a halo is described by the standard Navarro, Frenk, & White (NFW) density 
profile |5^, p oc (cr/i?200m)~^(l + cr/i?200m)~^, where -R200m is the halo virial radius 



which we define to contain a mean density of 200 times the mean density of the 
universe, the normalization is set by the fact that the mass profile must integrate to 
contain the total virial mass m within i?200m, and c is the halo concentration which 
sets the radius of the transition between the two power laws as ~ -R200m/c. 

The standard practice is to set the average concentration of a halo of mass m ac- 
cording to a phenomenological law derived from dissipationless A^-body simulations, 
such as 

c(m, z) = Co[m/m.,o]""(l + ^)"^ (2.7) 

where Co~10,a~0.1,/?~l ||5^, |5^, and m^o is the mass of a typical object 
collapsing at 2; = 0. We neglect the spread in halo concentrations at fixed mass 
because this is gives rise to a neghgible effect on the scales we consider, and on 
smaller scales this dispersion is degenerate with an overall shift in the concentration- 
mass relation [^, The shortcoming of Eq. ( p.7|) is that it describes halos in 
dissipationless A^-body simulations that neglect the physics of baryons. Fortunately, 
Ref. PT[ demonstrated that adopting a modified concentration-mass law within the 
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halo model suffices to model the convergence power spectra predicted by the baryonic 



simulations of Ref. to a level where biases in inferred dark energy parameters are 
less than 10% of their statistical uncertainties for forthcoming surveys. We use this 
as a starting point for our preliminary estimate of the influence of baryonic processes 
on tests of the consistency of general relativity through weak gravitational lensing. 

2.2 Parameterized Tests of the Consistency of General Relativity 

As we mentioned in the introductory section, the ideal scenario would be to test 
parameterized families of theories for modified gravity that make specific and unique 
predictions; however, this is difficult at present. Few if any viable alternatives to gen- 
eral relativity that contain a contemporary epoch of accelerated expansion have been 
identified. Moreover, weak lensing requires some treatment of the evolution of den- 
sity perturbations beyond the linear order of perturbation theory, and the nonlinear 
evolution of perturbations in models of modified gravity are not completely specified, 
nor have they been studied thoroughly. We consider two different parameterizations 
that have been proposed to test for the consistency of general relativity within forth- 
coming data. These do not represent complete, entirely self-consistent descriptions 
of any particular phenomenological alternative general relativity. We consider them 
as a pragmatic step toward testing the null hypothesis of gravity described by gen- 
eral relativity in the case of modest departures from the standard gravity. In both 
cases, we explore deviations to the linear growth of perturbations but assume that 
the relation between linear and nonlinear perturbations is unchanged. In practice, 
this means that we assume that halos of dark matter form with the abundances 
and other gross properties that they otherwise would have in a standard-gravity 
treatment and continue to employ the halo model to predict convergence spectra 
[Eq. on nonlinear scales (£ > 300). Modified gravity models require some en- 

vironmental dependence to the gravitational force law so that gravity may deviate 
from general relativity on large scales (low density) yet satisfy local constraints on 
deviations from general relativity (in the high density environment of our galaxy) 
[pH, Bl, IHl, B3, O . We do not consider any such modifications because there 



is no comprehensive treatment of the nonlinear evolution of perturbations in such 
theories and no phenomenological model akin to the halo model in the case of general 
relativity, though preliminary work in this direction has begun ^ . 



Our ffist test is to split the dark energy equation of state parameter into two 
distinct parameters. The ffist parameter, Wg, is used in all calculations of the growth 
of density perturbations while the second, Wd, is used in all calculations of the re- 
lationship between redshift and distance. The deflection of light is given by a sum 
of the Newtonian and curvature potentials and in the context of general relativity 
with zero anisotropic stress perturbation, these potentials are equal^. This equality 



^We are admittedly a bit cavalier here as there are several different notational conventions in use, 
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means that perturbation growth and the distance-redshift relation are both deter- 
mined by the evolution of the cosmic expansion rate H{z), so that Wg and Wd should 
be equal in general relativity and that an indication to the contrary may be a sign of 
modified gravity. Of course, any observational or theoretical systematic errors that 
drive the inferred values of Wg ^ Wd must be controlled and accounted for in order 
for this program to work. An alternative way to state our present aim is that we 
seek to assess the importance of the inability to predict the influence of the baryonic 
component of the universe on lensing power spectra for such consistency checks and 
explore a method to ensure the robustness of such a consistency check. In this first 
test, we consider our parameter set for dark energy to consist of the present dark 
energy density in units of the critical density Qbe, as well as the two equation of 
state parameters Wd and Wg. We take no priors on these parameters and set their 
fiducial values to f^oE = 0.76, Wd = uig = —1. 

In the second case, we explore the gravitational growth index parameter 7, in- 
troduced by hinder in Ref. |TT[ and explored in Refs. |T^, |T3|. hinder showed that 



a natural separation between the expansion history of the universe and the growth 
of perturbations could be achieved by taking the evolution of an overdensity 6 to be 
given by 

^ = ^uiay. (2.8) 
dm a 

Our notation is such that a = 1/(1 + 2;) is the cosmic scale factor, Qu with no 
argument of scale factor or redshift is the current density of non-relativistic matter 
in units of the critical density, ^Im{(i) = ^m(i~^ Hq / H [a)"^ with an explicit argument 
represents the evolution of the ratio of the matter density to the critical density, and 
Hq = H{z = 0) = lOOh km/s/Mpc is the present Hubble expansion rate. Eq. ( |2.^ ) 
holds for perturbations independent of scale or cosmic expansion history (aside from 
the dilution of f2M(a). In an enormous variety of dark energy models embedded 
within general relativity the growth index obeys 0, IT^, 



7 = 0.55 + 0.02[l + w(z = 1)], for w < -1 (2.9) 
= 0.55 + 0.05(1 + w(z = 1)], for w > -1, 

where w{z) is the dark energy equation of state evaluated at redshift z. Deviations 
from this relation therefore, may indicate non-standard gravity. This treatment 
is accurate to better than a percent in specific cases such as the self-accelerating 
braneworld gravity model of Ref. (so-called DGP gravity, which is now ruled 



see Refs. ^ [7^, ^ |7^; however, our contribution is not to explore the detailed behavior of 
the metric potentials in any specific theory (largely for lack of such options) , so we forego a detailed 
discussion. In fact, we have already assumed the equality of the Newtonian and curvature potentials 
in order to derive Eq. ( ^.l| ) for the convergence power spectrum in terms of the power spectrum 
of density fluctuations and which would otherwise be written in terms of the power spectra of the 
metric potentials directly. 
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out as an explanation of the cosmic acceleration, see Ref. but still serves as a 
useful example of the utility of this parameterization), which yields 7 ~ 0.69, as well 
as scalar-tensor and f{R) theories of gravity ||1T|, |T5[. In this case, we consider the 
expansion history to be dictated by dark energy with abundance ^de and a time- 
dependent effective equation of state w{a) = wo + Wa(l — cl) while perturbation 
growth is dictated by 7. We take no prior constraints on these parameters and 
perturb about the fiducial values f^DE = 0.76 (as in the first case), wq = —1, Wa_ = 0, 
and 7 = 0.55. 

2.3 Cosmological Parameters and Future Lensing Surveys 

In the previous section, we described the parameters we use to describe deviations 
from general relativistic gravity. To summarize, we take two test of the consis- 
tency of general relativity. In both cases, we assume dark energy with a present 
density of ^de = 0.76 in our fiducial model. In the first case, we split the dark 
energy equation of state parameter into two pieces, one governing the growth of den- 
sity inhomogeneities Wg, and one governing the distance redshift relation Wd- We 
take as a fiducial model Wg = = and assume no priors on either parame- 
ter. In the second case, we assume dark energy with an effective equation of state 
w{a) = Wq + Wa(l — «), with fiducial values Wq = —1 and Wa = 0. In this case we 
assume inhomogeneity growth parameterized by 7 according to Eq. ( |2.8| ), and take 
the fiducial value of 7 = 0.55 as in general relativistic gravity with a cosmological 
constant causing contemporary acceleration. 

Beyond the parameters that describe the dark energy /modified gravity sector, we 
consider four other cosmological parameters that influence lensing power spectra and 
may be degenerate with dark energy parameters. These parameters and their fiducial 
values are: the non-relativistic matter density um = ^uh'^ = 0.13, the baryon density 
ujb = flBh"^ = 0.0223, the amplitude of primordial curvature fluctuations = 
2.1 X 10~^ (in practice we actually vary In A^) at the pivot scale of kp = 0.05 Mpc~^, 
and the power-law index of the spectrum of primordial density perturbations rig = 
0.96. We adopt conservative priors on each of these additional parameters that are 
comparable to contemporary constraints on each of these parameters individually 
[0. To be speciflc, we take prior constraints of cr(ti;M) = 0.007, cr(ti;B) = 10~^, 
(T(lnA^) = 0.1, and o"(ns) = 0.04. Lastly, we allow three parameters that describe 
the effective concentrations of halos to be determined internally from the same data. 
These are the parameters cq = 13, a = 0.05, and /3 = 1 of the power-law halo 
mass-concentration relation in Eq. (|2.7|) . The flducial values we choose, Cq ~ 10 
and a ~ 0.1, differ from those found in N-body simulations of structure formation 
[discussion after Eq. (|2.7| )]. We choose these parameters to reflect the results of the 



hydrodynamic simulations of Ref. [p5 |. 



We study constraints that one would expect from three forthcoming galaxy imag- 
ing surveys. As mentioned in § |2.1|, we follow current convention by setting the intrin- 
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Figure 1 : The influence of the dark energy equation of state parameter on the relationship 
between angular diameter distance and redshift. The left panel shows the angular diameter 
distance as a function of redshift for models with tt^d 7^ — 1, in units of the angular diameter 
distance in our fiducial model with Wd = —1- The upper line shows -Da with Wd = —1.1 
and the lower line shows Da with Wd = —0.9. The right panel shows the net influence of 
this distance change on the lensing kernel of Eq. (|2.2| ) . The dashed line in this panel should 
be read against the right, vertical axis and represents the absolute weight W3{z)/ Hq. The 
solid lines show the particular example of the weight for the third tomographic bin ^3(2;) 
relative to its value in the Wd = case and should be read against the left vertical axis. 
The correspondence with iVd is as in the left panel. 

sic shape noise of source galaxies to be (7^) = 0.2 for all experiments and subsuming 
differences in the shape noise between surveys into differences in the effective number 
density of galaxies on the sky for each survey. We assume a redshift distribution of 
source galaxies as given in Eq. (p.3|). 

The most near-term survey that we explore is the dark energy survey (DES), 
which may begin operations in 2009 and have results as soon as 2011-2012^. For 
the DES, we take /sky = 0.12 and n = 15/arcmin^. We also consider an imaging 
survey that might be conducted as part of a future space-based mission such as 
the proposed Supernova Acceleration Probe (SNAP)^, which is a canonical example 
of a National Aeronautics and Space Administration, Beyond Einstein, Joint Dark 
Energy Mission probe^, though the specifics of the mission have yet to be decided and 
other competitors include the Advanced Dark Energy Physics Telescope (ADEPT) ^ 
and the Dark Energy Space Telescope (DESTINY)^. For a SNAP-like probe, we 

^http : / / www . darkenergysurvey . org 
^http : / / snap . Ibl . gov 

^http : / / universe . nasa . gov/pr ogram/ probes/ j dem 
^http : / / universe . nasa . gov/pr ogram/ probes/ adept 
^http : //www. noao . edu/noao/staff/lauer /destiny 



- 10 - 



1.04 
1.03 
1.02 

cri.oi 



Q 0.99 
0.98 
0.97 
0.96 

1 1.5 2 2.5 3 3.5 4 
1+z 

Figure 2: The influence of the dark energy equation of state parameter on the hnear 
growth of structure over an observationally-relevant range of redshifts for forthcoming 
imaging surveys. We plot the linear growth function for dark energy with Wg = —1.1 
(upper line) and Wg = —0.9 {lower line) as relative to the growth function in a cosmological 
constant model with Wg = —1.0. These relative shifts should be read against the left, 
vertical axis. In analogy with Fig. we plot the absolute growth function for the standard 
case of t(;g = — 1 as the dashed line, which should be read against the right, vertical axis. 
The dotted line shows the growth function in a flat cosmological model with 0,m = 1 and 
f^DE = (and should again be read against the right, vertical axis). Comparing this growth 
rate to the growth rate in the fiducial model with Ode = 0.76 illustrates the suppression 
of growth caused by the contemporary epoch of accelerated expansion. 




take /sky = 0.025 and n = 100/arcmin^. Lastly, we consider a future ground-based 
imaging survey as might be carried out by the Large Synoptic Survey Telescope 
(LSST)"^. We adopt /sky = 0.5 and n = 50/arcmin^ to describe the LSST survey. 
In all cases we consider only the statistical limitations of the surveys. The only 
systematic uncertainty we consider is the theoretical uncertainty associated with the 
inability to make precise predictions of the net influence of baryons on the lensing 
power spectra. 

3. RESULTS 

3.1 Consistency Checks 

We begin by illustrating the power of weak lensing surveys to test the consistency of 

^http : / /Isst . org 
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general relativity with forthcoming data. As we have mentioned above, the consis- 
tency checks work by comparing simultaneous constraints on cosmological distances 
and structure growth and determining whether they are consistent with a single, 
underlying, general relativistic model. Before moving on to parameter constraints 
we briefly discuss each of these effects individually. We will use the split dark en- 
ergy parameterization model (with parameters Wg and w^) to illustrate the influence 
of modified distances and modified growth. Those readers with significant experi- 
ence with dark energy phenomenology and proposed consistency checks of general 
relativity may like to proceed to the next section. 

The influence of dark energy on the cosmic distance scale is demonstrated in 
Figure |I[ The left panel of this figure shows the angular diameter distance as a func- 
tion of redshift in dark energy models with constant equations of state Wd = —1-1, 
Wd = —1.0, and = —0.9. The more negative the equation of state parameter, 
the more recent the prevalence of dark energy, and the more rapid is the current 
acceleration. Thus, for a fixed present Hubble expansion rate, lower values of Wd 
lead to larger angular diameter distances. The dependence of the angular diame- 
ter distance on the dark energy equation of state manifests in the lensing weight 
[Eq. ( p.2|) ]; this dependence implies that the lensing weight varies with Wd- The 
effects of this dependence are shown in the right panel of Fig. |I] for the particular 
case of the weight for the third tomographic bin containing W^lz), which contains 
sources with photometric redshifts in the range 1.2 < Zp < 1.8. Notice that the 
lens weight itself (dashed line) extends to 2; > 1.8. In our model, sources in this 
bin have true redshifts that extend beyond z = 1.8 because of the relatively large 
dispersion in calibrated photometric redshifts. For more negative Wd, the W^lz) is 
relatively increased and weighted relatively more toward low redshift. Heuristically, 
this can be understood by considering a hypothetical individual deflector. Relative 
to the deflector, the angle of deflection of a light ray is fixed by the deflector proper- 
ties and the apparent change of position of the source as seen by an observer grows 
with angular diameter distance. The varying distances to fixed redshift also affect 
the mapping between multipole and wavenumber in the matter power spectrum at 
fixed redshift. Increased distances mean smaller wavenumber at fixed redshift, and 
as the power spectrum of density fluctuations is a declining function of wavenumber 
on scales of interest, this results in slightly stronger lensing. 

Next, we demonstrate the dependence of the cosmic growth function on dark 
energy. Figure ^ shows the growth function D{z) in our fiducial cosmology (dashed 
line), normalized such that D{z) — (1 + z)~^ as z ^ 00. The dashed line in this 
figure shows the growth function in a flat cosmological model with no dark energy 
and flu = 1- Comparing these two growth functions shows the dramatic suppression 
in the growth of density perturbations during the epochs where dark energy makes a 
significant contribution to the cosmic energy budget. In addition, we show the growth 
functions in cosmological models with Wg = —1.1 and Wg = —0.9 relative to the 



- 12 - 



-0.4 




300 1000 3000 

I 



Figure 3: Partial derivatives of tomographic weak lensing power spectra with respect to 
cosmological parameters. We show derivatives of each of the power spectra P^^{t), with 
ij = 22, 23, 33 with respect to the three parameters w^, Wg, and In A^. The solid hnes show 
derivatives with respect to Wd- From bottom to top at left, the solid lines are OItiPk /dwd 
with ij = 33,23,22. The dashed lines represent derivative with respect to Wg. From 
bottom to top at left, the dashed lines show din /dwg with ij = 22,23,33. The dotted 
lines are derivatives with respect to the power spectrum normalization parameter InA^. 
From bottom to top at left, the dotted lines show —dlnP^/dlnA'j^ with ij = 22,23,33. 
The additional negative sign for the derivatives with respect to In is designed to reduce 
the dynamic range on the vertical axis. 

growth function in our fiducial model. The more positive the dark energy equation 
of state is, the earlier the acceleration begins and the earlier the growth of structure 
is quenched by the competing cosmic expansion. Conversely, more negative values 
of the dark energy equation of state lead to more recent expansion and, as a result, 
more aggregate growth of density perturbations since the epoch of recombination. 
These features are all represented in Fig. |[ 

The modifications to the growth of structure and cosmic distances manifest as 
changes in the observable power spectra. Figure |^ shows partial derivatives of three 
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of the power spectra, P^^, P^^, and P^^, as a function of parameters Wd, Wg, and 
In A^. The derivatives with respect to the dark energy parameters are all negative. 
To reduce the range of the vertical axis, we actually display — dP^j/dln A^, so that 
it will also lie in the negative vertical half plane. Notice that variations in the ob- 
servable power spectra due to changes in Wg and Wd have very different redshift and 
scale dependence. Consequently, we should not expect these parameters to exhibit 
a significant degeneracy and we would expect that shear power spectrum observa- 
tions could constrain both parameters independently. Notice also that changes in 
Wg induce relatively larger shifts in the observable power spectra than do shifts in 
Wd- In the absence of additional information, this might be taken as evidence that 
weak lensing constraints on Wg should be more stringent than such constraints on 
Wd- In practice, Wg is strongly degenerate with the power spectrum normalization 
parameter In A^. This can also be seen in Fig. |^. The dotted lines in Fig. ^ shows 
— dP^j/dln A^, and this quantity exhibits a similar redshift and scale dependence 
as the derivatives of the spectra with respect to Wg- Though not shown, Wg is also 
strongly degenerate with u^- The net result of these degeneracies is that Wd will 
be significantly more strongly constrained by weak lensing power spectrum observa- 



tions than Wg- This discussion is closely related to the eigenmode analysis of Ref. [16 



and serves as a qualitative demonstration that both Wg and Wd can be constrained 
by cosmic shear measurements, but the relative constraints on each parameter are 
sensitive to choices for external priors in the available additional parameter space. 

Having reviewed the influence of independent changes in the cosmological dis- 
tance scale and the rate of density fluctuation growth on observable weak lensing 
power spectra, we now turn to projections for the utility of general relativity consis- 
tency checks. We utilized the split dark energy equation of state parameterization 
for the previous illustrations, but we present results for both the Wd-Wg parameter- 
ization as well as the growth index (7) parameterization. Standard projections for 
the power of forthcoming consistency checks using weak lensing are summarized in 
Figure ^. Clearly such consistency checks can be quite powerful when definitive pre- 
dictions can be made for the convergence power spectra. In the particular case of 
LSST, projected la constraint on Wd is roughly a{wd) — 0.02, while the projected 
constraint on the gravitational growth index is cr{'~f) ~ 0.04, meaning that the LSST 
weak lensing program alone could rule out DGP gravity at the ~ 3.5(7 level. In the 
case of the Wd-Wg parameterization, it is the equality of these two parameters that 
serves as a null hypothesis, so it is useful to look at the marginalized constraint on 
the combination Wdm = Wd — Wg- Transforming to this variable, we find a la con- 
straint cr(w7difr) = 0.04 for LSST, a{wdis) = 0.11 for SNAP, and cr(«;difr) = 0.14 for 
DES. 

Unfortunately, these standard projections assume perfect knowledge of nonlinear 
structure growth and ignore the influence of poorly-understood baryonic processes on 
the predictions of lensing spectra. Uncertainties due to baryonic processes and the 



-14- 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

"■-1.3-1.2-1.1 -1 -0.9-0.8-0.7 -1.3-1.2-1.1 -1 -0.9-0.8-0.7 

Figure 4: The capability of forthcoming surveys to perform consistency checks of general 
relativity. The left panel shows projected la constraint contours in the Wd-Wg plane for 
the LSST (innermost contour), DES (outermost contour), and SNAP weak lensing (middle 
contour) experiments. The diagonal line in this panel delineates Wg = Wd. The constraints 
on the parameter lUdiff = — Wg from each experiment are given in the lower, right 
portion of the panel. The right panel shows constraints in the gravitational growth index 
parameterization in the wq-^ plane. The three contours represent the LSST, DES, and 
SNAP weak lensing experiments as in the left panel. 

process of galaxy formation are significant and will be the subject of the remainder 
of this section, beginning with the following section where we demonstrate that such 
effects can significantly bias estimators of dark energy parameters and diminish the 
utility of consistency checks if not accounted for. 

3.2 Biases 

The poorly-constrained baryonic processes are detrimental to any consistency check 
of general relativity because, if unaccounted for, they may induce biases in inferred 
cosmological parameters. In particular, if the estimators of dark energy parameters 
that we have introduced are significantly biased, it would be possible to conclude 
erroneously that data are inconsistent with general relativity. The relevant quantities 
to examine in order to assess the importance of baryonic effects are the biases in 
units of the statistical uncertainties of those parameters. If the biases are small 
compared to the statistical uncertainties, it will be unlikely to rule out the true 
model based on these biases, but as biases become comparable to or larger than 
statistical uncertainties it becomes increasingly likely that the true model may be 
ruled out based on the data. 

In practice it is difficult to estimate what biases may be realized, because the true 
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Figure 5: Biases in the estimators for the spUt dark energy parameters Wg and u^d that 
may be reahzed if baryonic processes are ignored. The upper, left panel shows the the effect 
of modified halo structure on convergence power spectra. Each of the lines that increases 
with multipole represents the relative change in a convergence power spectrum of sources in 
our third tomographic bin (1.2 < Zp < 1.8) P^^{i), in models with modified halo structure 
relative to that of the standard case. We represent the standard case using a halo model 
where halo concentrations are given by Eq. ( p.7| ) with cq = 10, a = 0.1, and P = 1.0. 
We represent models with modified halo structure by taking cq = 11, 12, 13, 14, and, 15 
from bottom to top. The shaded bands show the statistical errors on P^^(^) expected 
from forthcoming SNAP, DES, and LSST surveys from top to bottom at left. The other 
three panels show the biases in estimators w^^ (solid) and Wg (dashed), in units of the 
statistical uncertainties in these parameters, as a function of the maximum multipole used 
in parameter estimation. These should be read against the left vertical axes. Each panel 
shows forecasts for a specific survey. The most biased cases correspond to cq = 15. The 
decreasing functions of ima,x in each panel show the statistical uncertainty in Wd and Wg as 
a function of ^max relative to the error if all information to £ma.x = 3000 were used. These 
lines should be read against the right vertical axes. 
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power spectra must be known in order to perform such a calculation (and, of course, 
if the true spectra were calculable there would be no bias). Fig. |^ gives an estimate 
of biases that may reasonably be realized for a variety of assumptions regarding the 
true convergence power spectra. We constructed these bias estimates as follows. Fol- 
lowing the simulation analysis of Refs. we assume that the primary influence 
of baryons can be accounted for in a halo model by adopting a non-standard halo 
concentration relation. It has been demonstrated that incorporating this additional 
freedom into the halo model enables accurate modeling of the results of non-radiative 
and dissipational hydrodynamic cosmological simulations and reduces expected bi- 
ases in inferred dark energy parameters to acceptable levels . The standard result 
for the concentration relation from gravity-only simulations is Eq. (|2.7|) with roughly 



Co = 10, a = 0.1, and (3 = 1.0 |59|. We take this model to be our treatment of the 



results of A^-body simulations. The halos in the simulations of Ref. |]35[ exhibit a 
~ 40% enhancement in their concentrations in dissipational simulations relative to 
A^-body simulations, so a choice of cq ~ 14 is an acceptable model for these simu- 
lation results. Unfortunately, our understanding of the evolution of baryons is poor 
and no simulation campaign can produce definitive results (it is this lack of definitive- 
ness that is the primary problem). In particular, this boost is likely an over-estimate 
of any true shift due to galaxy formation because contemporary simulations exhibit 
baryonic cooling and star formation that is more efficient than allowed observation- 



ally fTj, |7|, 0. As a consequence, baryonic processes must be treated in a 

way that accommodates a wide range of predictions. We illustrate reasonable levels 
of bias due to galaxy formation effects by computing the biases induced in dark en- 
ergy parameter estimators when the convergence spectra are modeled with a fixed 
halo concentration- mass relation with Cq = 10, but where the "true" value of this 
normalization ranges among cq = 11, 12, 13, 14, and, 15. The lower end of this range 
is only a small shift relative to the standard concentration relation and is realized 
even in non-radiative hydrodynamic simulations where baryons cannot radiate their 
energy and form galaxies. The upper end of this range is slightly larger than the 
Ref. simulation results and is near the maximum allowable by current lensing 
constraints at low redshift {z ~ 0.2) We note in passing that there are as 

yet weak observational indications of concentrations higher than those predicted by 
dissipationless simulations [0, |8| . 

Fig. ^ shows biases in the inferred Wg and Wd induced by treating convergence 
power spectra that would be well described by an enhanced concentration mass 
relation, where Cq can be any of Cq = 11, 12, 13, 14, or 15, with a concentration-mass 
relation fixed to have Cq = 10. The objective is to provide reasonable estimates for 
biases that may be realized by not treating baryonic processes appropriately. If the 
Ref. simulation results were a good representation of the "true" convergence 
power spectra, the biases would be near the top of the range in Fig. |^, roughly 
corresponding to the fourth most strongly biased model in each panel. 
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To see why it is useful to consider biases in dark energy parameter estimators as 
a function of scale may require some elaboration. The influence of baryons is greater 
on smaller scales, so a simple way to eliminate the uncertainty caused by baryons is 
to excise small-scale information from parameter estimation. A simplistic way to do 
this is to choose a maximum multipole £max for the analysis and to disregard spectra 
beyond that maximum (Ref. explores more sophisticated means to excise small- 
scale information, but they yield similar results in this context |^|). In Fig. ^, we 
show biases as a function of £max in order to illustrate the reduced biases realized as 
a function of minimum scale (maximum-multipole). 

Fig. ^ illustrates several points of interest. Most previous lensing analysis have 
assumed that one could utihze information to at least £max ~ 10^ effectively. Fig. ^ 
shows that if one considers such small-scale information, estimators of dark energy 
equation of state parameters Wg and will both be significantly biased relative to 
their true values. Notice also that they are generally biased in an opposing sense. This 
is of paramount importance for consistency checks where the goal is not necessarily 
to constrain either parameter individually, but to test the hypothesis that Wg and 
Wd are the same. For example, if Wg and were always biased in the same way due 
to an unaccounted for systematic then it would not be possible to produce a reliable 
constraint on either parameter, but in this contrived example it would still be possible 
to test for their equality without accounting for the systematic. The biases in Fig. ^ 
are generally large compared to statistical uncertainties and suggest that consistency 
checks of general relativity can be compromised by uncertainty in galaxy formation 
physics and may lead to rejection of the hypothesis of general relativity even in 
models where this is the true theory of gravity. The point where biases become 
acceptable {S^ <^ cr^) varies from one experiment to the next. One could argue that 
for the LSST survey these biases are not acceptable for any £max ^ 300. 

As a rough criterion. Fig. ^ demonstrates that reducing the biases induced by 
the shifts in halo structure to acceptable levels requires taking a maximum multipole 
of no more than a few hundred. Of course, this reduction in information comes 
at a cost. The monotonically decreasing lines in Fig. |^ show the marginalized la 
statistical uncertainties in Wd and function of the maximum multipole £max 

in units of the uncertainty when all information to an i^nax = 3000 is used. Notice 
that these uncertainties should be read relative to the right vertical axes in each 
panel. Overlaying the dependence of parameter uncertainties on this plot serves to 
demonstrate the cost of excising small-scale information. For LSST, the cost is large, 
a factor of ~ 3 — 4 in the uncertainties of both Wd and Wg, because of the exquisite 
precision with which LSST could measure convergence spectra. This greatly reduces 
the effectiveness of an LSST-like data set to test the consistency of general relativity. 
DES is least affected by the uncertainty in galaxy formation because it surveys a 
comparably low number density of galaxies and makes the least use of small-scale 
information of any of the three experiments. SNAP is intermediate between the two 
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Figure 6: Estimates of biases in the estimator of the 7 parameter of the gravitational 
growth index formahsm proposed in Ref . 1 11 ] . The treatment of baryonic modifications to 
the convergence power spectrum is as in Fig. ^ and the upper, left panel is identical to the 
upper, left panel of Fig. |5|. We repeat it here for convenience. The remaining panels show 
biases projected for the LSST, SNAP, and DES weak lensing surveys as indicated. Value 
of the biases should be read relative to the left vertical axes. The biases are given in units 
of the projected statistical uncertainty in 7 and the most (least) biased case corresponds 
to a model with cq = 15 (cq = 11). The intermediate lines correspond to the intermediate 
values of cq with bias increasing with cq. The monotonically decreasing lines in each 
panel represent the la statistical uncertainties in 7 as a function of maximum multipole 
^max- These are shown in units of the uncertainties in gamma achieved by considering all 
information to ^max = 3000 in order to illustrate the relative degradation in constraints 
caused by excising small-scale information. 



with constraints degraded by roughly a factor of ~ 2.5. 

Figure |^ shows the biased estimator for the gravitational growth index introduced 
in Ref. 111]]. Our methods for computing biases and the presentation are analogous 



to that of Fig. ^ Of course, in this case the only parameter of importance is the 
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gravitational growth index 7. Qualitatively, Fig. ^ shows results that are similar 
to that of Fig. ||. In particular, biases can be considerable compared to statistical 
uncertainties. Moreover, excising small-scale information comes at great cost in the 
statistical uncertainties of the dark energy parameters, with relative degradation of 
~ 2 — 4 depending upon the choice of imax and experimental parameters. 

3.3 Calibration 

In the previous section, we demonstrated the significance of the biases that would 
be induced in tests of general relativity if baryonic processes were to be ignored. We 
also showed the significant cost of excising information from the small scales where 
baryonic effects are most important in terms of degraded parameter uncertainties and 
thus degraded tests of the theory of gravity. What remains is to study whether such 
tests can still be performed utilizing high-multipole information if baryonic effects are 
treated in a parameterized way and calibrated simultaneously with the cosmological 
parameters. Though there is still significant work to be done, preliminary indications 
are that baryonic effects can indeed be treated by modifying the predictions of dissi- 



pationless A^-body simulations ||35|, |3^. In fact, the necessary modification is to treat 
the relation between halo mass and halo concentration as free, and we have already 
utilized such modifications to estimate parameter biases in the previous section. 

The next step is to determine the degree to which treating both halo structure 
and cosmology as uncertain, in order to eliminate the biases of the previous section, 
and fitting both to forthcoming data will degrade our ability to test general rela- 
tivity. If the degradation is significantly less than that incurred by simply excising 
small-scale information, then an appropriate strategy to adopt in the run-up to the 
next generation of imaging surveys would be to spend significant effort understand- 
ing the phenomenology of structure growth in the nonlinear regime in a variety of 
different effective models for the evolution of the baryonic sector. It might then be 
possible to "self-calibrate" the influence of baryons, simultaneously yielding insight 
into the law of gravity on large scales and the small-scale physics of galaxy forma- 
tion. Alternatively, if the degradation is comparable to the degradation incurred by 
excising small-scale information, such a program would likely be futile and it would 
be wiser to excise small-scale information until such time as definitive predictions for 
the influence of baryonic processes on lensing power spectra can be made. 

There is already reason to be optimistic that some self-calibration can be achieved. 
Ref. and Ref. |29| demonstrated the robustness of weak lensing to modest mul- 



tiplicative and additive errors in shear measurements. In addition, Ref. have 
already demonstrated robustness of weak lensing constraints on dark energy to some 
uncertainty in the power spectrum at high wavenumbers. A large part of potential to 
carry out self-calibration programs within weak lensing data can be gleaned from the 



eigenmode analysis of Ref. ||T6[. The robustness of weak lensing constraints on dark 



energy is due in part to the fact that weak lensing provides exquisite constraints on 
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distance measures that are relatively insensitive to dark energy and can be brought 
to bear on other systematics without degrading dark energy constraints. 

In this section, we present the results of a self-calibration exercise where we 
treat cosmological parameters as before and introduce additional parameter freedom 
to describe the influence of baryons. To be specific, we treat the halo concentrations 
according to Eq. {\i.7\), but allow the parameters Cq, a, and P to be free. It is already 
possible to offer a guess at the promise of such a self-calibration exercise. Fig. ^ 
shows the derivatives of three of our convergence power spectra P^^, P^^, and P^^ 
with respect to the parameters w^, Wg, and In A^. Compare this to the shifts in P^^ 
depicted in the upper, left panels of Fig. ^. Modifying the concentration relation 
causes significant increases in power at small scales (high multipoles). As a result, 
there is a stronger scale dependence associated with changes in concentration param- 
eters than there is associated with changes to Wd, Wg, 7, or A^. This is encouraging 
because it suggests that if the scale dependence can be adequately modeled, the 
parameters describing baryonic processes can be extracted independently without 
significantly degrading consistency checks of general relativity. 

We show the results of the self-calibration exercise in Figure ^ In this panel 
we show constraints in the case of self calibration of halo structure cr**'^'^~'^^'(£max), 
in units of the parameter uncertainties that would be quoted in the standard case 
of perfect knowledge of halo structure, cr^^^'^'^^^'^iimax) (shown in Fig. ^, Fig. |^, and 
Fig. In particular, we show how the parameter degradation varies with £max 
because it is both instructive and it illustrates the balance between excising data 
and self-calibrating. Figure |^ shows results for both the w^-Wg parameterization as 
well as the gravitational growth index (7) parameterization. 

It may seem odd that the functions in Fig. |^ do not decrease monotonically 
with increasing £max- These relations are not monotonic functions of £max because 
they represent the ratio of the error realized in the self-calibration calculation com- 
pared to the errors computed in the limit of perfect knowledge of the influence of 
baryons. From Fig. ^ and Fig. ^ the constraints on all parameters are rapidly de- 
creasing functions of ^max- This is sensible; as more information is added constraints 
should improve and not degrade. Likewise, in the case of self-calibration the absolute 
constraints on each parameter decrease rapidly with increasing ^max- What is not 
monotonic is the relative degradation of these constraints as a function of ^max- 

The scale dependence of the parameter degradation is not surprising. From 
the power spectra in the upper, left panel of Fig. |^ and the derivatives in Fig. || 
it is evident that each of the parameters of interest induces a scale-dependence on 
the observed spectra, so some scales are more effective at constraining dark energy 
parameters and calibrating halo structure than others. Though the details depend 
upon the experiment, and in particular the statistical weight an experiment places 
on a particular multipole, there are sensible general trends depicted in Fig. |^ In 
the limit of low-£, halo structure is unimportant and results in the self-calibration 
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Figure 7: Results of self-calibrating halo structure simultaneously with tests for modified 
gravity on large scales. As with Fig. ^, we present results as a function of the maximum 
multipole considered in the parameter extraction analysis £max- The vertical axis is the error 
on the stated model parameter in the self-calibration case in units of cr^"-"-'^"-^'^ : the error 
on parameter X calculated by standard parameter forecasts that appear in the literature, 
where halo structure is assumed to be known perfectly. Each panel contains three lines for 
each set of experimental parameters we explore. The solid lines correspond to an LSST-like 
experiment, the dashed lines correspond to a SNAP-like experiment, and the dotted lines 
correspond to a DES-like experiment. Each panel focuses on a different parameter. The 
upper, left panel shows the upper, right panel shows Wg, the bottom, left panel shows 
the difference Wdis = Wd — Wg, and the bottom, right panel shows the gravitational growth 
index 7. Note that Wd and Wg are the parameters of our "dark energy split" model and 
are constrained simultaneously, and that the relevant constraint for tests of gravity is the 
constraint on Wdis- The 7 parameter is from an entirely distinct parameterization in which 
dark energy is independently marginalized over. The additional lines that separate this 
panel from the others are intended to reinforce this distinction. 

calculation approach the standard results. This is least evident for LSST because 
LSST has sufficient /sky to make precise power spectrum measurements even for 
i ~ 300. We do not extend these plots to i < 300 because at this value of £ 
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Figure 8: Constraint contours in the W([-Wg plane after accounting for and marginalizing 
over the uncertainty in halo structure. The filled, grey contours correspond to the la 
contours in this plane from the LSST-like (innermost contour), SNAP-like (middle contour), 
and DES-like (outermost blue contour) experiments we consider. The thick, solid lines 
overlaying the contours represent the confidence contours in the limit of perfect knowledge 
of nonlinear structure formation as shown in the left panel of Fig. |^. The cross shows the 
fiducial model and the diagonal line running from the lower left to the upper right delineates 
W(]^=Wg. As an explicit demonstration of the effect of excising small-scale information to 
eliminate potential bias, the outermost, blue contour shows the la marginalized constraint 
in the Wd-Wg plane from LSST if no information from multipoles greater than £max = 300 
are utilized. 

parameter biases are already at or near acceptable levels and so there is no need 
to relegate ourselves to such low multipoles in any data analysis. With increasing 
multipole, the halo structure becomes ever more important and is mildly degenerate 
with dark energy. Further increasing imax beyond ~ several hundred introduces the 
strong scale dependence that halo structure imparts on the convergence spectra and 
this breaks a degeneracy between dark energy and halo structure. For this reason, 
the degradation decreases rapidly from i 

max 

500 to i 

max 

2000. 

The most important aspect of Fig. |^ is that the degradations are relatively small 
if we utilize all of the information to £max ^10'^. Degradations are scale dependent, 
but consider the level of degradation at imax = 3000 for definiteness. After calibrating 
halo structure, the Wd constraint is degraded by ~ 34% for LSST and roughly ~ 18% 
for DES and SNAP. The Wg constraint is degraded by ~ 18% for LSST and less than 
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10% for LSST and SNAP. The relevant combimation Wdis = Wd — Wg is degraded 
by a maximum of ~ 36% for LSST and only about ~ 13% for SNAP and DES. To 
appreciate the utility of self- calibration in this instance, we should compare these 
numbers to the information loss associated with excising small scale information in 
Fig. 1^. For LSST, excising high-£ information to eliminate biases requires taking 
^max ~ 300 with a corresponding ~ 300% increase in the constraints on and Wg. 
So long as baryonic effects can be modeled in this way, the more sensible strategy is 
clearly to use all available information and calibrate the physics of galaxy formation 
as encoded in effective halo concentrations. The rightmost panel of Fig. ^ shows a 
similar result for the gravitational growth index 7. In particular, is degraded by 
only ~ 20% after self calibration as compared to ~ 200% as would be required by 
the excision of small scales (see Fig. |^). 

For completeness, we give a revised set of constraint contours in the w^-Wg plane 
in Fig. ^. After calibration, the contours are slightly expanded, but their orienta- 
tion in this plane is only slightly changed. By comparison, eliminating small-scale 
information corresponds to a dramatic loss of constraining power both because the 
area of the constraint contour in the w^-Wg plane grows and because its orientation 
changes such that the most degenerate combination becomes more nearly perpindic- 
ular to the line = Wg. The overall conclusion remains that self calibration is the 
appropriate strategy to adopt. 

4. Discussion 

We have extended recent studies of the influence of galaxy formation on weak lensing 
observables (in particular, Ref. [^) to address cross-checks of the consistency of gen- 
eral relativity with weak lensing. Weak lensing is a key ingredient in such programs 
because weak lensing observables are sensitive to both the cosmological distance 
scale and the evolution of potential inhomogeneities, whereas the Type la super- 
novae distance-redshift relation and baryonic acoustic oscillation measurements are 
sensitive only to geometry (though it may be possible to exploit the lensing of super- 
ovae while performing the distance-redshift test, see Refs. ^ ^] or to combine 
distance measures with cluster counts). Proposed methods to utilize weak lensing 
to constrain dark energy and to place limits on deviations from general relativity on 
cosmological scales assume that correlations on scales as small as a few arcminutes 
{i ~ 10^) can be brought to bear on the cause of cosmic acceleration. Such small- 
scale information is important to these programs. Relegating consideration to larger 
scales can significantly degrade the constraining power of forthcoming experiments 
(see Fig. |^ and Fig. ||) by nearly a factor of ~ 3 decrease in the effectiveness of tests 
for deviations from general relativity. 

One of the anticipated challenges of using this small-scale information is that it 
is difficult to predict observables on these scales because the physics that governs 
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the evolution of the baryonic component of the Universe is poorly understood and 
affects these observables in important ways ||3^, |33|, ^ In Fig. | and Fig. H, 
we show that using such data to test general relativity may introduce potentially 
large biases in the inferred values of cosmological parameters that can lead to ruling 
out general relativity when it is in fact the true theory of gravitation. In these 
figures, we give a range of estimates of potential biases that may be realized. These 
results are motivated by the specific results of the simulations of [^. The biases 
can be many times the statistical uncertainties in these parameters, which is clearly 
a serious problem. We reiterate here that we give a range of biases that may be 
realized precisely because it is difficult to predict the precise influences of baryons. 
We have used simulation results as guidance and explored a range of concentration 
relations that are consistent with current observations (e.g., Ref. p2[). We have 
not included information from additional observables that can be brought to bear 
on this problem. For example, supernova la and baryon acoustic oscillation data 
can constrain the cosmic distance-redshift relation and be combined with lensing 
to produce stronger consistency checks. However, these additional constraints on 
cosmological distances would only make the discord between structure growth and 
distance measures more egregious and drive larger biases, making the problem more 
severe. 

Faced with a large systematic error, a common practice is to degrade the experi- 
ment in some way so as to derive robust results from the data and reduce sensitivity to 
the systematic error. Baryonic influences are scale dependent, so one way to degrade 
the experiment and eliminate these biases is to disregard small-scale information by 
restricting the analysis to £max ~ a few x 10^. However, in Fig. ^ we have shown that 
this comes at the cost of nearly a factor of ~ 3 decrease in the effectiveness of tests 
for deviations from general relativity (more sophisticated analyses, such as those in 
Ref. |^5|], yield similar results). 

As an alternative to excising small scale information, simulations suggest that 
we can understand the form of the baryonic influence and correct for it using a 
parameterized model. To the degree that this form is not exactly degenerate with 
parameterizations of dark energy or deviations from general relativity, this enables 
the utilization of some of the information on small scales. We addressed these biases 
by including 3 baryonic parameters, and Cq, that model the effective concen- 
tration of dark matter halos as specified in Eq. |2.7| . This choice of parametrization 



is motivated by Refs. fSg, wherein the authors demonstrated that modifying 



halo concentrations allows one to model account for the net infiuences of baryonic 
processes on scales relevant to weak lensing parameter estimation. In particular, 
models in which the A^-body results are corrected by a modified halo concentra- 
tion relation track the convergence spectra from galaxy formation simulations out to 
^ ~ 5000 accurately enough to reduce the biases in inferred dark energy parameters 
to levels below the statistical uncertainty of forthcoming experiments. Within this 



- 25 - 



framework, we studied the ability of forthcoming imaging surveys to self-cahbrate 
the values of the baryonic parameters and simultaneously use this information to 
check the consistency of general relativity. 

In Fig. we have illustrated the effectiveness of our method by using the Fisher 
Matrix formalism to estimate the constraints on 7 and Wg — Wd that will be obtained 
from future photometric surveys. From this figure, it is evident that the constraints 
on the consistency of general relativity degrade by just ~ 30% relative to the limit of 
perfect knowledge of nonlinear evolution of the gravitational potential. In Fig. |^ the 
constraints on the consistency of general relativity obtained through our methods 
are much tighter than those from analyses excising small scale information. In light 
of our findings, we conclude that self-calibration of some of the specific influences 
of baryonic physics may be an appealing alternative to disregarding the information 
contained on scales smaller than £ ~ 300. 

Our outlook for self-calibration of uncertain baryonic physics is quite positive; 
however, there are several important caveats to our study and a number of addi- 
tional studies that must be undertaken in order to bring such a program to fruition. 
First, this correction for the physics of galaxy formation has only been applied to 
the simulations of Ref. [^. As the physics of baryons is highly uncertain, a more 
comprehensive exploration of viable alternative models is needed in order to validate 
such a correction term. In fact, it is not reasonable to suppose that the influences of 



baryons should be strictly confined to the internal structures of halos (see Ref. [39 



for changes in the halo mass function) and other effects may prove important. Nev- 
ertheless, our study indicates that such an exploration is a fruitful pursuit. 

Second, any such tests have important limitations. The self-calibration program 
that we have explored should guard against the possibility of ruling out general 
relativity when it is the correct theory of gravity. In our particular study, we have 
not explored complete, self-consistent phenomenologies for modified gravity. Rather, 
we have only explored proposals that are already in the literature in which additional 
parameters are introduced which have known, fixed values in general relativity and 
explored the ability of forthcoming surveys to constrain these parameters. In the 
event of a detection of a deviation from general relativity such tests provide limited 
information about the character of the correct, alternative theory. In large part, this 
is due to the fact that modified gravity models must have some scale dependence so 
that they can simultaneously drive accelerated cosmic expansion yet satisfy small- 
scale bounds on the theory of gravity [^, ^ |6^, |6^, If a modified gravity 
theory, including nonlinear evolution, were fully specified and explored and if galaxy 
formation processes were more completely understood, it would be possible to draw 
more specific conclusions from forthcoming imaging survey data. 

Lastly, we note that we operate under the assumption that dissipationless A^- 
body simulations will effectively calibrate the properties of the lensing field absent 
baryonic effects. In our halo model approach, this is tantamount to a calibration of 
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dark matter halo abundance and halo clustering that is so accurate as to introduce 
negligible uncertainties. This is a common premise to all such studies of dark energy 
constraints as it is known that the density field is not yet calibrated with sufficient 
accuracy to analyze forthcoming data 0]. Our outlook is that it is plausible 

that a brute force simulation campaign can achieve sufficient accuracy in this regard, 
though it is certainly possible that reaching this goal may be complicated by several 
factors that have not yet been thoroughly studied. For example, the phenomenolo- 
gies we explore here are relatively simple, as they treat modified gravity only as a 
modification to the growth function. Ultimately, it will be of interest to study more 
specific and more complete phenomenological models of modified gravity that may 
lead to large corrections to the density field on scales of several Mpc. Additionally, 
even in the absence of corrections necessitated by modified gravity it is possible that 
a program of iV-body simulations coupled with a calibration of baryonic effects may 
encounter unforeseen obstacles. In either case, the assumption that an adequate 
suite of A^-body simulations will address the density field (in our parameterization 
halo bias and halo abundance) with sufficient precision so as to render errors in the 
predicted density field negligible in the dissipationless case may not be valid. In such 
cases, it would be necessary to allow for adequate model freedom to account for this 
uncertainty. In the language of our paper, this might amount to internal calibration 
of halo abundance and halo bias (presumably with significant prior information) as 
well. Such additional freedom would lead to further degradation in dark energy 
and/or modified gravity parameters. We have chosen to limit the scope of our work 
along these lines and not studied this more complicated case. 



5. Summary 

We have explored the ability of forthcoming weak lensing surveys to perform a con- 
sistency check on general relativity as the correct theory of gravity in light of recent 
uncertainties in predicting lensing observables on scales smaller than ~ 10 arcmin 
due to the poorly-understood evolution of the baryonic component of the Universe. 



After a pedagogical demonstration of the utility of such consistency checks in § 
we present the following results. 



1. Conducting a study of the consistency of general relativity by analyzing weak 
lensing observables from a forthcoming imaging survey (such as DES, JDEM/SNAP, 
or LSST) out to scales as small as £ ~ 10^ could potentially lead to inferred cos- 
mological parameters that are biased by many times the statistical uncertainty 
of such surveys if baryonic effects on the nonlinear evolution of the gravita- 
tional potential are ignored. The problem is severe in that they could lead to 
rejection of the null hypothesis of general relativistic gravity, even when it is 
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the true theory. Specific results depend upon the assumed survey properties 
and the importance of baryonic processes. 

2. Disregarding small-scale information beyond £ ~ a few x 10^ significantly re- 
duces these biases because of the strong scale-dependence of the infiuence of 
baryonic physics on weak lensing observables. However, excising this infor- 
mation degrades the constraining power of forthcoming imaging surveys by a 
factor of ~ 3. 

3. As an alternative to excising small-scale information, we explore a program of 
self- calibrating a phenomenological model of baryonic processes that is moti- 
vated by recent numerical simulations. With such a model small-scale informa- 
tion out to £ ~ 10^ can be exploited, and biases can be reduced to levels below 
those of the statistical uncertainty of forthcoming experiments. As shown in 
Fig. ^ and Fig. |^, this method provides substantially tighter constraints on 
the consistency of general relativity than analyses that disregard small-scale 
information. 

4. The program to calibrate the baryonic correction internally allows small-scale 
information to be used to constrain gravity and guards against ruling out gen- 
eral relativity when it is the true theory. In the case of internal calibration, 
the limits provided by general relativity consistency checks are degraded by 
only ~ 30% relative to the limits that would be achieved if nonlinear structure 
formation were perfectly understood. Moreover, interesting limits on galaxy 
formation models are produced "byproduct." 

These results should be carefully considered in preparation for the weak lensing 
surveys that will be conducted in the next decade. At the very least, it should be 
clear that issues regarding the baryonic influences on lensing observables will need to 
be dealt with, but that uncertainty in this regime may not necessarily have dramatic 
consequences on the ability of imaging surveys to illuminate the cause of cosmic 
acceleration. It seems likely that a program of internal calibration of baryonic physics 
will allow lensing surveys to achieve stringent limits on gravity and dark energy. In 
order to ensure that this is the comprehensive theoretical program will need 

to be undertaken to better understand viable alternatives for the net influences of 
baryons on observable power spectra. 
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